On the isoscalar-isovector splitting of pygmy dipole structures 
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The electric dipole response of ^'"'Ce is investigated using the fully consistent relativistic quasi- 
particle random phase approximation. By analyzing the isospin structure of the El response, it 
is shown that the low-energy (pygmy) strength separates into two segments with different isospin 
character. The more pronounced pygmy structure at lower energy is composed of predominantly 
isoscalar states with surface-peaked transition densities. At somewhat higher energy the calculated 
El strength is primarily of isovector character, as expected for the low-energy tail of the giant 
dipole resonance. The results are in qualitative agreement with those obtained in recent (7,7') and 
(a, a'7) experiments, and provide a simple explanation for the splitting of low-energy El strength 
into two groups of states with different isospin structure and radial dependence of the corresponding 
transition densities. 
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Exotic modes of collective excitation represent unique 
structure phenomena in nuclear many-body systems with 
a pronounced asymmetry in the number of protons and 
neutrons A representative example is the pygmy 

dipole resonance (PDR), a low-energy mode that corre- 
sponds to a resonant oscillation of the weakly-bound neu- 
tron skin against the isospin saturated proton-neutron 
core. Even though high-resolution data on low-energy 
dipole states have been obtained in recent (7, 7') exper- 
iments d H, i, S, 0, H, the details of the underly- 
ing nature and dynamics of the PDR remain largely un- 
known. Recent experiments with radioactive ion beams 
have extended the studies of low-lying dipole response 
towards exotic nuclear systems with pronounced neutron 
excess [^, The occurrence of PDR has also been 

investigated using a variety of theoretical approaches 
and models fH. More recent studies have made use of 
the quasiparticle RPA plus phonon coupling [ll|, the 
quasiparticle-phonon model [l^l , the relativistic RPA ISj, 
14 1 and QRPA [lEj, as well as the relativistic quasiparti- 
cle time blocking approximation 1^, 17 1. 



A recent {a,a'^) study of electric-dipole excitations 
below the particle threshold in the semi-magic nucleus 
^^°Ce 0, has shown that the PDR structure separates 
into two parts. By comparing the observed low- lying 
strength to results from photon scattering experiments on 
the same target Q, it was noticed that the lower part of 
the PDR is excited both in [a, a'^) and (7, 7') scattering, 
whereas the higher-energy component seems not to be 
excited by a-particles. 

In this letter we analyze the isospin properties and 
transition densities of low-lying dipole excitations in 
^^°Ce. Similar aspects of the PDR will be considered 
as in Ref. , but from a theoretical point of view. The 
goal is to explore in simple terms, i.e. by analyzing the 
isotopic and spatial dependence of transition densities for 



dipole excitations, the splitting of the PDR structure into 
two groups of states, excited either by both the nuclear 
part of the a-nucleon interaction and photons, or only 
photons, respectively. 

The present analysis uses the fully self-consistent 
relativistic quasiparticle random phase approximation 
(RQRPA) based on the Relativistic Hartree-Bogoliubov 
model (RH B) Il8| . Details of the formalism can be found 
in Refs. [H, [if^In the RHB-I-RQRPA model the effec- 
tive interactions are implemented in a fully consistent 
way. In the particle-hole channel effective Lagrangians 
with density-dependent meson-nucleon couplings are em- 
ployed llil , and pairing correlations are described by the 
pairing part of the finite-range Gogny interaction [20| . 
Both in the ph and pp channels, the same interactions are 
used in the RHB equations that determine the canonical 
quasiparticle basis, and in the matrix equations of the 
RQRPA. 

The full set of RQRPA equations is solved by diago- 
nalization. The result are excitation energies E\ and the 
corresponding forward- and backward-going amplitudes, 
AT^ and Y^, respectively, that are used to evaluate the 
transition strength: 
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where fi and /i' denote canonical single-nucleon states, 
and M, V are the corresponding occupation factors. Dis- 
crete spectra are averaged with a Lorentzian distribution 
of arbitrary width (r=1.5 MeV in the present calcula- 
tion), and the electric El response is calculated for the 



2 




E[MeV] 



FIG. 1: The total RHB+RQRPA El transition strength for 
^*°Ce (sohd). The dotted, dashed, and dot-dashed curves 
connect states with predominant (at least 70%, 80% and 90%, 
respectively) isoscalar (IS) components, identified by analyz- 
ing the corresponding proton and neutron transition densities 
over the radial interval [0, 8] fm. The vertical arrow denotes 
the one-neutron separation energy Sn ~ 9.2 MeV. 



isovector dipole operator: 
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Particularly relevant for the present study are the neu- 
tron and proton radial transition densities 5pn,p{^:)^ which 
describe nucleon density fluctuations induced by an ex- 



ternal field (2l|. For instance, by analyzing if 5yOn(r) and 



5pp(j) for a particular state appear to be in phase over 
some extended region within the nuclear volume, states 
with predominant isoscalar structure can be identified in 
the RQRPA dipole strength distribution (Eq. (P)). 

The solid curve in Fig. 1 represents the total 
RHB+RQRPA El transition strength for ^'^'Ce, calcu- 
lated with the DD-ME2 [23] parameterization of the ef- 
fective Langrangian. In addition to the pronounced col- 
lective isovector giant dipole resonance (GDR) peaked 
at « 15 MeV, we notice an enhancement of El strength 
in the low-energy region below neutron threshold, indi- 
cated by the vertical arrow. This structure could be at- 
tributed to a PDR. The dotted, dashed, and dot-dashed 
curves illustrate the isospin structure of the El strength 
distribution. These curves connect states with predomi- 
nant (at least 70%, 80% and 90%, respectively) isoscalar 
components, identified by analyzing the corresponding 
proton and neutron transition densities over the radial 
interval [0, 8] fm. When for a particular state the proton 
and neutron transition densities are found to be in phase 
over more than 70% of this range of the radial coordi- 
nate, the state is denoted IS 70%, and analogously for the 
states IS 80% and IS 90%. For the low-energy excitations 



FIG. 2: Same as in Fig. 1, but the dotted, dashed, and dot- 
dashed curves are determined from the radial dependence of 
the neutron and proton transition densities in the surface re- 
gion: 4.25 < r < 8 fm. 



around the pronounced peak at 8.4 MeV the proton and 
neutron transition densities arc in phase over more than 
80% of the radial interval [0, 8] fm, i.e. these states are 
predominantly isoscalar. Some isoscalar states are also 
found at excitation energy E>10 MeV but in that in- 
terval, obviously, isovector collective motion dominates. 
The analysis of the isospin structure of the El transi- 
tion strength in Fig. 1 indicates that most of the low- 
energy strength is of isoscalar character, and therefore 
sensitive to isoscalar probes such as a-particles. States at 
higher energy, however, would probably not be excited in 
(a, a'7) experiments because of their predominant isovec- 
tor character. Photon scattering, on the other hand, 
should excite dipole states over a broader low-energy in- 
terval, but (7, 7') experiments have been restricted to the 
region below the neutron separation energy. Indeed, re- 
sults of RQRPA calculations [3, [l^ indicate that part of 
the PDR strength may be missing in (7, 7') experiments, 
i.e., these calculations predict PDR strength around and 
also above the neutron threshold. Some progress in re- 
covering part of the unresolved El strength has already 
been reported , but future experimental studies across 
the neutron threshold may be necessary to account for 
the missing low-lying El strength. 

The photon field interacts with a nucleus as a whole, 
whereas a-particles interact mainly with the nuclear sur- 
face, inducing isoscalar transitions with surface-peaked 
transition densities. It is, therefore, interesting to ana- 
lyze in more detail the radial structure of the transition 
densities for the dipole states shown in Fig. 1. When the 
identification of states with predominant isoscalar com- 
ponents is based on the behavior of the corresponding 
neutron and proton transition densities in the radial in- 
terval [0,4.25] fm, i.e., in the volume of the nucleus, one 
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FIG. 3: The RQRPA neutron and proton transition densities 
for the peaks at 8.4 MeV and 15.1 MeV excitation energy in 
i«Ce. 



finds that these states are distributed over a wide en- 
ergy region, but not below E « 10 MeV. Fig. 2, on the 
other hand, displays the same analysis of the El tran- 
sition strength, but the predominantly isoscalar states 
are identified by considering the in-phase or out-of-phase 
behavior of the corresponding neutron and proton tran- 
sition densities in the surface region 4.25 < r < 8 fm. 
In this case transitions with predominant isoscalar char- 
acter are located mainly in the low-energy region below 
10 MeV. One can conclude that the low-energy isoscalar 
structure peaked at 8.4 MeV corresponds mainly to sur- 
face excitations. Low-lying isoscalar states with surface- 
peaked transition densities are precisely those expected 
to be induced in a-scattering experiments. 

In Fig. 3 we plot the neutron and proton transition 
densities for the peaks at 8.4 MeV and 15.1 MeV. For 
the main peak at 15.1 MeV the transition densities dis- 
play a radial dependence that is characteristic for the 
isovector GDR: the proton and neutron densities oscil- 
late with opposite phases. The dynamics of the state at 
8.4 MeV is different: the proton and neutron transition 
densities are in phase in the nuclear interior and there 
is almost no contribution from the protons in the sur- 
face region. This low-energy state does not belong to 
statistical El excitations sitting on the tail of the GDR. 
Fig. 4 explores in more detail the structure of the transi- 
tion matrix elements for the PDR state at 8.4 MeV. The 
matrix elements (i.e. each term in the sum over fi and 
/i' in Eq. (1), representing the contribution of the cor- 
responding 2qp configuration to the transition strength) 
are plotted as a function of the unperturbed energy of the 
corresponding 2qp configurations. Proton and neutron 
matrix elements are displayed separately. First, we note 
that the state at 8.4 MeV is characterized by a coherent 
superposition (matrix elements of the same sign) of more 
than ten neutron 2qp configurations. The largest ma- 
trix elements correspond to the following neutron tran- 
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FIG. 4: The RQRPA neutron and proton transition matrix 
elements for the state at 8.4 MeV in ^'"^Ce, as a function of the 
unperturbed energy of the corresponding 2qp configurations. 



sitions (in order of increasing 2qp unperturbed energy): 

3si/2 —>■ 3p3/2, 2d3/2 3p3/2, 3si/2 —>■ 3pi/2, 2^3/2 ^ 
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from the anal- 



l/iii/2 25g/2, I57/2 2/5/2. 
ysis of the QRPA amplitudes (^2gp)^-(^29p)^ it follows 
that the total contribution of proton 2qp excitations to 
the state at 8.4 MeV is only 9%, well below the ratio 
Z/N expected for a GDR state, whereas 91% of the total 
amplitude corresponds to neutron transitions. 

The relativistic QRPA accurately reproduces experi- 
mental excitation energies of the GDR, whereas the cal- 
culated PDR are only in qualitative agreement with avail- 
able data. In the particular example of ^'^''Ce, the pre- 
dicted PDR structure is calculated more than 2 MeV 
above the pygmy structure seen in (7,7') and {a,a'j) 
experiments [2]. The principal reason for this discrep- 
ancy lies in the fact that relativistic mean-field effective 
interactions, e.g. DD-ME2, have relatively low effective 
nucleon masses, typically m* < 0.7 rriN, where ttin de- 
notes the free nucleon mass. A larger effective mass, i.e. a 
higher density of states around the Fermi surface, would 
lower the excitation energy of states with predominant 
isoscalar components, e.g. the PDR. Unfortunately, it 
is not simple to increase the effective nucleon mass in 
relativistic mean-field models, because this quantity is 
strongly constrained by the Dirac mass, which represents 
a measure of the strength of the spin-orbit single-nucleon 
potential [23^, ^2^ . In addition, our study is limited to a 
simple QRPA and does not include effects of coupling to 
low-energy surface phonons. This coupling will also en- 
hance the effective nucleon mass. The analysis of PDR 
can be extended by introducing the coupling of two- 
quasiparticle excitations to collective vibrations, e.g. us- 
ing the quasiparticle phonon model [l3 | , or the relativis- 
tic quasiparticle time blocking approximation flij. In 
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FIG. 5: The total RHB+RQRPA El transition strength for 
^^''Ce in comparison with calculation using experimental neu- 
tron single-particle energies. 



Ref. [16| it has been shown that the effect of two-phonon 
admixtures is a small shift of PDR states (< 1 MeV) to 
lower excitation energies but, although the PDR calcu- 
lated in the extended space contains sizable two-phonon 
admixtures, it basically retains a one-phonon character 
and the pygmy dynamics is not modified by the coupling 
to low-energy surface vibrations. 

To illustrate how sensitive PDR is to the effective nu- 
cleon mass, in Fig. 5 we display the fully self-consistent 
RQRPA El strength calculated with the DD-ME2 inter- 
action, in comparison to the El strength obtained with 
the same effective interaction, but with DD-ME2 neu- 
tron single-particle energies below and above the iV = 82 
shell gap replaced with experimental values. The empir- 
ical = 82 gap is much smaller than the one predicted 
by DD-ME2, i.e. the empirical effective neutron mass 
is larger. By using experimental neutron single-particle 
energies one takes into account, at least on the level of 
single-nucleon energies, the effect of coupling to collective 
vibrations but, of course, the calculation is no longer self- 
consistent. Nevertheless, we notice that the enhancement 
of the effective neutron nucleon mass has a pronounced 
effect especially on the PDR structure. The El strength 
below 10 MeV is fragmented into two structures, with the 
lower one centered at « 5 MeV and higher one peaked 
at « 8 MeV, in much better agreement with data 0. 
An analysis of the isospin content of these structures, 
based on proton and neutron transition densities, shows 
that isoscalar components are more pronounced in the 
lower structure at w 5 MeV. We have also verified that 
the same pattern is calculated for ^^^Sn, when DD-ME2 
single nucleon energies are replaced by empirical values. 

In conclusion, the fully consistent relativistic QRPA 
has been used to calculate the electric dipole response of 
^^°Ce. An analysis of the isospin structure of El strength. 



based on the radial dependence of the microscopic pro- 
ton and neutron transition densities, has shown that the 
low-energy strength separates into two segments with 
qualitatively different isospin character. The more pro- 
nounced PDR structure at lower energy is composed of 
predominantly isoscalar states with surface-peaked tran- 
sition densities. At energies above the main PDR peak 
the El strength is primarily of isovector character, with 
non-surface peaked transition densities, as expected for 
the low-energy tail of the giant dipole resonance. The 
results are in qualitative agreement with those obtained 
in recent (7,7') and {a,a'^) experiments, and provide 
a simple explanation for the splitting of low-energy El 
strength into two groups of states, excited either by both 
the nuclear part of the a-nucleon interaction and pho- 
tons, or only photons, respectively. 
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